Metal-insulator transition in EuO 
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It is shown that the spectacular metal-insulator transition in Eu-rich EuO can be simulated within 
an extended Kondo lattice model. The different orders of magnitude of the jump in resistivity in 
dependence on the concentration of oxygen vacancies as well as the low-temperature resistance mini- 
mum in high-resistivity samples are reproduced quantitatively. The huge colossal magnetoresistance 
(CMR) is calculated and discussed. 



I. INTRODUCTION 

Stoichiometric (pure) EuO is a ferromagnetic semicon- 
ductor. In the sixties and seventies it became famous for 
the redshift of its optical absorption edge below the Curie 
temperature Tc (69 K)i. Eu-rich EuO, whose Eu rich- 
ness is realized by oxygen vacancies, behaves like a metal 
at low temperatures and shows a tremendous metal- 
insulator transition near the Curie temperatur o 2 ' 3 ' 4 . A 
jump in resistivity of 13 orders of magnitude was mea- 
sured by Penney et al^ but, since leakage currents limited 
the measurement, it was probably even greater. It is a 
remarkable feature that the size of the jump in resistivity 
varies greatly amongst different experimental samples 3 . 
Another intriguing feature is a low-temperature resis- 
tance minimum in high-resistivity samples. Further- 
more, Eu-rich EuO has a huge colossal magnetoresistance 
(CMR)&. 

Recently interest in EuO was renewed by measure- 
ments of the resistance and of the spin-split conduction 
band of a EuO film by Steeneken et alA The authors 
emphasize the extraordinary properties of Eu-rich EuO 
and their importance for basic research in the field of 
spintronics due to the almost perfect spin polarization of 
the conduction electrons at low temperatures. 

The first physical picture (He model) for the metal- 
insulator transition in Eu-rich EuO was developed 
by Oliver et al. 2 . They postulated a temperature- 
independent electronic trap (impurity) level created by 
the oxygen vacancies. It is below the conduction band for 
high temperatures, where Eu-rich EuO is an insulator. 
On decreasing the temperature below Tc, the spin-up 
conduction band is shifted towards lower energies (red- 
shift) and therefore crosses the impurity level at some 
point. The electrons of the trap level are emptied into the 
conduction band giving rise to EuO becoming a metal. 
Later the model was refined by introducing a Coulomb 
repulsion between the two electrons at an oxygen vacancy 
and ascribing a temperature-dependent shift to the trap 
levels^. However, the shift of the impurity levels is un- 
physical as will be explained below. There is a second 
physical picture (the bound magnetic polaron (BMP)), 
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which was introduced by Torrance et ali They assumed 
an exchange interaction of the surplus electrons of each 
vacancy site with the Eu spins. Theoretical attempts to 
quantify the metal-insulator transition of Eu-rich EuO 
were made by Laks et al£, Mauger— and SteenekeniS. 
The conductivity in all three cases was obtained by the 
simple Drude formula a = n c efj, (n c - conduction-electron 
density, e-elementary charge, /^-mobility) . 

In this paper we formulate an extended Kondo-lattice 
Hamiltonian as a microscopic model for the metal- 
insulator transition in Eu-rich EuO. It contains qualita- 
tively the idea of the He model through the assumption 
of impurity trap levels, and it contains the physics of the 
BMP since an exchange between the impurity-electron 
spins and the Eu spins is included. However, our formu- 
lation goes beyond the former two approaches because 
of the detailed modelling of the impurity orbital energies 
and their hybridization with the conduction band. The 
disorder of the oxygen vacancies is treated within an ex- 
tension of the coherent-potential approximation (CPA). 
Our calcu lation s should be more accurate than the ones 
in Refs. Isi lillOl since we are using a more complex con- 
ductivity formulai 1 * 1 ^, which is derived directly from the 
fundamental current-current correlation Kubo formula. 



II. THEORY 

In order to model the situation in Eu-rich EuO, first 
one has to consider a conduction band eg. The spins of 
the conduction electrons <Xj are coupled to the Eu spins 
Si via a coupling constant J. These terms represent the 
Kondo lattice model. Then, there are impurity levels e p , 
a Coulomb repulsion U in case of two electrons at an 
oxygen vacancy, and an exchange between the impurity- 
electron spins a? and the Eu spins with coupling constant 
J p . Furthermore, there should be a hybridization V be- 
tween impurity electrons (Pj a ) and conduction electrons 
(ct CT ). The superexchange between the Eu spins is not 
modeled explicitly. Instead we are considering (S z ) as a 
parameter, whose temperature dependence is given by a 
Brillouin function, which is fulfilled to very high accuracy 
in Eu-rich EuQi 3 -. The Hamiltonian which captures all 



2 



the features discussed is the following: 

ka » 

3 a 3 3 

The sums over i and j mean a sum over all lattice sites 
and a sum over the randomly distributed oxygen- vacancy 
(impurity) sites, respectively. As the actual value of J p 
is not known, J p = J is assumed for simplicity. 
The current operator is given by 



3 = lim j^[H,p{q}] 
q^o \q\ 



(2) 



with the density operator 
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P(<1) = 



k — qa ka / j 

ka 3° 



p} aPj a + 0(q<). (3) 



Eqs. ©-© yield a surprisingly simple result for the cur- 
rent operator in a cubic system: 



V 9k x ka ko 



(4) 



Eq. Q has to be put into the current-current correlation 
Kubo formula for the conductivity^. If the self-energy 
is local, which is exact in infinite dimensions, e.g. if ap- 
plying dynamical- mean field theory (DMFT), there is a 
substantial simplification to the Kubo formula. All ver- 
tex corrections to the current-current correlation function 
vanish and one is left with an expression that only con- 
tains the one-particle spectral densityii*i£. We apply the 
local approximation to the self-energy and use the same 
conductivity formula for the three-dimensional case: 



,/£l-/'(E)a A .jE) 2 (^ 



ka -co 



(5) 



where V is the volume, f'(E) is the derivative of the 
Fermi function and A?(E) is the spectral density of the 
conduction electrons. Eq. (JSJ) can be transformed^ into 

oo oo 

a = £E/ dE (-f'( E )) J dx<j>{x,E)v{x) (6) 
<f>(x,E)=A j ; a (E) 2 



v(x) = - / dE'E' Po (E') 



a is the lattice constant (= 5.1 A in EuO), and po(E) is 
the free conduction band density of states. 



To get the spectral density A^ (T (E), we introduce ap- 
propriate selfenergies: 

«[<%,. - J E^" *<]-: C L» = ^a(E) «(*„; 4» (7) 

i 

3 i 
= Y<l{E)((p ra -p), a )) . (8) 

((...;...}} stands for the retarded Green's function. 
Eqs. and JSJ are so far exact relations. 

The main approximation of our approach consists in 
making independent ansatzes for the conduction-electron 
and impurity selfenergies. This is justified if the effect of 
the hybridization, which to the lowest order is propor- 
tional to V 2 , is small. In our calculations it turns out 
that the hybridization V itself has to be small to give 
a reasonable order of magnitude of the metal-insulator 
transition. Therefore, it is reasonable to neglect in the 
selfenergies effects of the order of V 2 . 

The conduction-electron self-energy is taken from an 
interpolating ansatz 17 for the conduction electron part 
of the Hamiltonian (Yl e k n ka ~~ ^Yl^i ' 

ka i 

2 4 1 - b a Go(E + p - ±z a J{S z )) 

(9) 

a a = S(S + 1) - z a {S z ){z a {S z ) + 1) , K = \j 



where z a — 5 a ^ — <5 CT j and Gq(E) is the free conduc- 
tion electron Green's function. The ansatz fulfills all 
known limiting cases for n c — > (ferromagnetic satura- 
tion, atomic limit, second-order perturbation theory in 
J and high-energy-expansion up to the fourth moment). 
It is therefore especially appropriate in the present case 
of very small conduction-electron densities (10~ 3 - 1CP 4 
per unit cell). 

For the impurity self-energy Y%(E) we take the self- 
energy of the atomic limit of the correlated Kondo- 
lattice modeli& (impurity part of the Hamiltonian Q) 

(E e P n i<r + f/ E r 4r n il ~ J pY.Sj ■ Oji)- If u is lar g e 

3<? j j 

enough in comparison with J, this leads to a 4-peak struc- 
ture in the impurity quasiparticle density of states with 
the highest peak at £4 = e p + U + ^S. The peak posi- 
tions are fixed, i.e. do not depend on the magnetization 
(S z ). Therefore, the above-mentioned picture by Oliver 
et al£, which asigns a temperature-dependent shift to 
the impurity levels, is wrong. However, the weights of 
the peaks depend on (S z ) as well as on the impurity oc- 
cupation. For any parameters only 3 peaks have a finite 
weight. The peak at £4 will act as the impurity level 
which donates electrons to the conduction band. 

We have made ansatzes for the selfenergies Tj p a (E) and 
Yi a {E) which are correct for all known limiting cases of 
n c — > and up to the first order in the hybridization V. 
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To treat the randomness of the oxygen vacancies we 
use the coherent-potential approximation (CPA), which 
originally was invented to deal with a non-interacting 
alloy 19 . Here wc arc dealing with strongly interacting 
electrons in a two-component alloy. One component cor- 
responds to the impurity sites A, the other to the non- 
impurity sites B. First, because of our effective medium 
approach and the locality of the impurity self-energy, lo- 
cal effective one-particle energies for the impurity sites 
are retained. The second problem is how to treat the 
non-impurity sites in a consistent way. To solve this, for 
non-impurity sites we also introduce one-particle levels 
X) e p,i n ia- an d hybridization terms with the conduction 

a 

band V J2(Pia c i<? + AoVio) ■ The one-particle energies e Pii 

a 

must diverge to infinity at non-impurity sites to assure 
that these levels are never really occupied. With these 
virtual levels and a non-random hybridization on the lat- 
tice the usual CPA equation can be formulated: 







H + Y? m<T {E)-Y.<f\E) 



i=A.B 



l-G p a {E) (e p , m - m + S^OE) 



£CPA ( £)) 

(10) 



e p m = A 
oo m = B 



YP a m = A 
m = B 



with ca = d and cb = 1 — d the concentrations of impu- 
rity and non-impurity sites, respectively. is the local 
impurity Green's function and S^ PA the CPA self-energy. 
The CPA is known to be correct up to the first order in 
the impurity concentration The Green's function of 
the conduction electrons is given by 



1 



E - (e(k) - fi) - £ „ s jjPA (g) - Za(E) 



(11) 



From Eq.^Jthe spectral density A^ a (E) can be obtained. 
Since Aj: (T (E) 2 enters the conductivity formula (JSJ), our 
calculations of the resisitivity are correct for all limit- 
ing cases of n c — > and up to the second order in the 
hybridization V and the impurity concentration d. 

The parameters for the calculations are as follows. The 
width of the conduction band W (« lOeV) can be taken 
from the absorption spectrum of Ref. or from band 
structure calculations 21 :. We have assumed a semiellip- 
tical shape of the conduction band, which approximates 
well the actual situation at the lower band edg e2i. J de- 
termines the shift of the conduction band (jS), which is 
about 0.3 eV (half band-splitting) at low temperatures^. 
For Eu spins of S = | this implies a coupling constant 
J = 0.17eV. The position of £4 is in general not known. 
Only for the highest-resistivity curves an activation en- 
ergy of 0.3 eV was observed^. If U is assumed to be 
fixed at leV, the impurity energy e p remains the param- 
eter to be adjusted. The hybridization V, which is also 



not known, should not be too large in order not to destroy 
the effect of the metal-insulator transition. We have cho- 
sen V = O.OlcV. Of all parameters, W and J are taken 
from the experiment. The precise values of U and V are 
not important apart from the fact that the hybridiza- 
tion should be small. The decisive free parameter is the 
impurity level £4 since it determines the position of the 
chemical potential. For the vacancy (impurity) concen- 
tration d medium values are of the order of 0.1 %££ Two 
electrons are assigned to each oxygen vacancy. 



III. RESULTS AND DISCUSSION 
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FIG. 1: Comparison of measured (diamonds, samples 66-6, 
95-BA-3 and 49-4 from Ref: 3 ) and calculated resistivity (solid 
lines) in dependence on temperature. The theoretical curves 
are calculated for three different impurity concentrations d. 
€ p = — 6.44eV (t4 = — 5.14eV, lower band edge at — 5.30eV at 
T = 0). For other parameters see the text. 

In Fig. ^ a comparison of calculated and measured re- 
sistivity curves is shown. The measured curves represent 
three examples of moderate resistivity from Ref. y- The 
only varying parameter of the three theoretical curves is 
the vacancy concentration d. e p (and with it 64) was fixed 
so as to yield the best overall agreement with the exper- 
imental curves. A variation of e p (£4) has to be taken 
into account in order to explain the resistance minimum 
in high- resistivity samples (see below). Each curve in 
Fig. shows a huge metal-insulator transition of sev- 
eral orders of magnitude near the magnetic transition 
temperature. The resistivity in the whole temperature 
range gets higher and the jump in resistivity at T c big- 
ger, the fewer oxygen vacancies (electron donors) there 



4 



are. Note the different orders of magnitude (7, 6 and 
5) of the jump in resistivity of the three experimental 
samples, which together with the absolute values of the 
resistivity are reproduced fairly well by the calculated 
curves. This demonstrates that the variation in the re- 
sistivity behaviour can be explained by different impurity 
concentrations. Incidentally, there is almost a 100% spin 
polarization of the conduction electrons at low tempera- 
tures as reported by Steeneken et al£, which is important 
for possible applications in spintronics. 



0.25 
0.2 

0.15 
0.1 

0.05 







- ! 1 0.04 


* — 






0.04 




s 
i 


- ! i 0.02 


/Hi 
" / : i- 






0.02 






- !t=iok] 

- i 10 

" ! 

- 1 
I 


/ :l: 




- |t=60K 







iiL u- 


-0.1 0.1 




" i 
- i 
i 
j 


-0.1 


0.1 


- j 

i. ,, ! . / 


/ 

', i 




- i 
_ i 
i 




/ ✓ 

f t 
t , 


1 1 


2 -1 ( 


) 1 2 - 


2 -1 





l : 



0.25 r- 

0.2 - 
0.15- T=/0K 

0.1 - 
0.05 - 



0.04 
0.02 




-0.1 0.1 



-1 



FIG. 2: Conduction-electron densities of states 

(Pt( — )>P|( )) an d impurity densities of states 

(p P t(- ■ -),Ppl(---)) for e p = -6.44eV, d = 0.05% and 
three different temperatures. For other parameters see the 
text. The insets show the densities of states near the Fermi 
level. 

A qualitative explanation of the metal-insulator tran- 
sition is possible in terms of quasiparticle densities of 
states along the lines of the He model (with the differ- 
ence that there are three instead of one impurity level, so 
the number of electrons which can be emptied into the 
conduction band is not so big). As shown in Fig. [3 for 
low temperatures the uppermost impurity level £4 lies 
within the spin-up part of the conduction band, which 
is down-shifted due to the exchange with the Eu spins. 
Therefore, a part of the impurity electrons is emptied 
into the conduction band. The system is a metal. On in- 
creasing the temperature the conduction band is shifted 
upwards. The number of conduction electrons decreases 
and the resistivity rises. Above the Curie temperature 
the conduction band is well separated from the upper- 
most impurity level, so only thermally excited electrons 
contribute to the conductivity, like in an n-doped semi- 
conductor. The resistivity is high and decreases with 
increasing temperature. Hence, the system is an insula- 
tor. 

For a quantitative analysis of the metal-insulator tran- 



sition we rely on the conductivity formula Eq. JJJ. The 
energy integrals in that formula yield a highly non-linear 
dependence of the conductivity on the number of conduc- 
tion electrons n c , which goes beyond the simple Drudc 
formula a = n c efi applied by other authors. 

For high-resistivity samples there is a characteristic 
low-temperature minimum in the resistivity^. Oliver et 
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FIG. 3: Comparison of measured (squares, sample 34-2-30 
from Ref£) and calculated resistivity (solid line) for a high- 
resistivity sample. Impurity concentration d — 0.01%, e p — 
— 6.61eV (t4 = — 5.31eV, lower band edge at — 5.30eV at T = 
0). For other parameters see the text. 

al£ explained it by activation from additional impurities 
like Gd or La. However, a more direct explanation is 
possible within our model without assuming additional 
impurities. The high values of the resistivity indicate 
that the impurity peak at £4 must be rather distant from 
the conduction band at Tc- On lowering the tempera- 
ture, the resistivity first decreases since the spin-up con- 
duction band moves down towards the impurity level £4. 
More electrons are thermally excited from the impurity 
level into the conduction band. If the position of £4 is 
low enough never to cross the conduction band, then on 
lowering the temperature further, the resistivity increases 
again near the ferromagnetic saturation {S z } « S. When 
the conduction band stays almost fixed with respect to 
the impurity level, fewer electrons are excited into the 
conduction band with decreasing temperature. To test 
this explanation, we have calculated the temperature- 
dependent resistivity for e p — — 6.61eV (£4 = — 5.31eV, 
lower band edge at — 5.30eV at T = 0) and impurity 
concentration d — 0.01%. In Fig. [31 the calculated and a 
measured* curve for a high-resistivity sample are shown. 
The theoretical curve has a minimum at T — 20K which 
fits quite well to the experimental one. Moreover, there is 
a quantitative agreement between the high-temperature 
tails of the experimental and the calculated curves. For 
the temperature range in between no measured points 
are available but the run of the theoretical curve seems 
credible although it takes on very high values. (Penney 
et al.'s above-mentioned resistivity measurements^ were 
limited to values of 10 n f2cm but interpolating their data 
allows values of up to 10 16 f2cm.) 

The calculated dependence of the resistivity on a mag- 
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netic field B is shown in the left graph of Fig. 0] It 
looks qualitatively similar to the measured dependence 
in Fig. 3 of Ref. Oi A characteristic shift of the resistivity 
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FIG. 4: Left: Calculated resistivity for different magnetic 
fields B; right: negative magnetoresistance. B = 0T( — ), 

0.05T (•••), 0.1T ( ), 0.5T (— •), IT ( ), 2T ( ). 

d = 0.1%; for other parameters see the text. 

maximum from Tq towards higher temperatures can be 
observed. This is due to the shift of the conduction 
band depending on the magnetization. Applying a mag- 
netic field has a similar effect as lowering the tempera- 
ture. At the same time the value of the maximum de- 
creases significantly since the metal-insulator transition 
is smeared out because of the higher temperature of the 
maximum. Therefore there is a huge negative magnetore- 
sistance MR = p ^ B p(~B^ as shown in the right graph of 



Fig- IH If one normalizes the magnetoresistance with p(0) 
instead of p(B), the magnetoresistance will be practically 
1 in a wide temperature range. This can be compared 
with a value of only 0.8 for the colossal magnetoresistance 
(CMR) of the manganites like Lai-^Ca^MnO^. 



IV. CONCLUSIONS 

In summary, a model has been presented which, in con- 
nection with a conductivity formula that is derived from 
the fundamental Kubo formula in the local-self-energy 
approximation, allows to accurately reproduce the spec- 
tacular metal-insulator transition and the huge CMR in 
Eu-rich EuO. The oxygen vacancy concentration is the 
decisive parameter for the big variations in resistivity 
behaviour. It may influence the position of the upper- 
most impurity level, which is important in order to ex- 
plain the resistance minimum in high-resistivity samples. 
The precise dependence of the resistivity on both the 
impurity concentration and the position of the impurity 
level should be a subject of thorough future experimental 
investigation, especially if aiming at applications in the 
field of spintronics. 
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